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ABSTRACT 

We present UBVRI photometry for the old open cluster NGC 6253. From comparison 
of the observed colour-magnitude diagrams with simulations based on stellar evolu- 
tionary models we derive in a self consistent way reddening, distance, and age of the 
cluster: E(B-V)= 0.23-0.32, (m-M) = 10.9±0.1, metallicity roughly double than so- 
lar, and age ~ 3 Gyr. The cluster has a binary sequence, discernible even through 
the field contamination, suggesting that about 1/3 of the cluster members belong to 
binary, or multiple, systems. 

Key words: Hertzsprung-Russcl (HR) diagram - open clusters and associations: 
individual: NGC 6253 - Age - Metallicity 



1 INTRODUCTION 



Open clusters are a key probe of the chemical evolution and 
dynamics of the Galaxy (e.g., Janes 1979, Panagia and Tosi 
1981, Friel and Janes 1993, Tosi 1996) and the interactions 
between thin and thick disks (e.g., Sandage 1988), since they 
provide more reliable information on the average stellar ages 
and radial velocities at different galactic radii (e.g., Janes 
and Phelps 1995) than field stars. They are among the very 
few galactic objects for which meaningful distances can be 
derived over a large range, which makes them an essential 
tool to constrain galactic evolution theories. However, to 
appropriately exploit the potentialities offered by open clus- 
ters, homogeneous and high quality observational material, 
as well as accurate data treatments, are mandatory to avoid 
misleading effects (e.g, Carraro and Chiosi 1994, Friel 1995). 

We have undertaken a project aimed to secure photo- 
metric observations of the required accuracy for a well cho- 
sen sample of open clusters mapping a grid of ages, metallic- 
ities and galactocentric distances. The goal is to expand the 
number of clusters for which good quality photometric data 
are available, to constrain current theories on the dynamics, 
structure and chemical evolution of the Galaxy. Relevant 
cluster parameters are usually obtained by fitting theoret- 
ical isochrones to the observed color-magnitude diagrams 
(CMDs). We follow a different approach, namely the com- 
parison of observed CMDs with synthetic ones generated by 
numerical codes based on stellar evolutionary tracks (Tosi et 
al. 1991), which we have found to be quite more powerful in 



Based on observations made at ESO telescopes, La Silla, Chile 



investigations both of galactic clusters (Bonifazi et al. 1990) 
and of nearby irregular galaxies (Marconi et al. 1995) . 

In this paper we examine NGC 6253, an old open cluster 
so far scantly studied, located toward the galactic center 
(Qi950 = 16 fe 55.1 ro ,5 195 o = -52°38'; l u = 336°, 6 n = -6°). 
Our results on Collinder 261 (age > 7 x 10 9 yrs and solar 
metallicity) have been presented by Gozzoli et al. (1996), 
while studies of NGC 2506 and NGC 6603 are in progress. 

In section 2 we describe the observations and data 
analysis; in Section 3 we present the derived CMDs 
involving U,B,V,R,I filters and discuss the problem of 
back/foreground contamination and the presence of a con- 
spicuous population of binary stars. Section 4 will be de- 
voted to the comparison with synthetic CMDs and the 
derivation of metallicity, age, distance and reddening, and 
the findings will be discussed in Section 5. 



2 OBSERVATION AND DATA REDUCTIONS 

We have observed NGC 6253 at the Danish 1.54m telescope 
located in La Silla, Chile. The CCD mounted in direct imag- 
ing was a Tek 1000x1000, with a pixel scale of 0.377 arc- 
sec/pixel and a total field of 6.8x6.8 arcmin. We observed 
two slightly overlapping fields, one centered at the cluster co- 
ordinates, and the other just North of it, intended as exter- 
nal field. Standard CCD fields (PG1323-086, PG1633+099, 
SA110, Mark A, T Phe; Landolt 1992) containing a total of 
18 standards both blue and red were also acquired, as was 
the usual set of bias, dark exposures, and sky flats. Expo- 
sures were taken in the Johnson-Cousins U,B,V,R,I filters, 
both very short (~ 10 seconds) to avoid saturation at the 



© 0000 RAS 



2 Bragaglia et al. 



Figure 1. Map of the observed region of NGC 6253; North is up and East left, and the field of view is about 7x7 arcmin. This map 
derives from our photometry, so it may be incomplete (e.g., there is a very bright field star at ~ 950, 270, not visible here). Also indicated 
are the four stars whose pixel and absolute coordinates are given in Table 2. 



Table 1. Journal of observations. Exposure times are given is seconds. 



Field 


Date 


U 


B 


V 


R 


I 


Central 


Jul 11-18, 1993 


900,900,120,40 


900,900,30 


600,600,300, 


600,120,120, 


600,600,300,180, 










60,10 


30,10,5 


60,30,30,5 


North 


Jul 11-18, 1993 




60 


30 




30 



bright end, and longer to reach deep on the main sequence. 
Table 1 gives a journal of the observations and the exposure 
times for all niters, while Figure |l| shows the central field 
here analysed. Table 2 lists the X, Y coordinates in pixels 
of four stars, and their equatorial coordinates, for possible 
conversion. The celestial coordinates were found identifying 
the stars in the Digitized Sky Survey^] images distributed on 



CD-ROM. Precision is about 1 arcsec, both in right ascen- 
sion and declination. 

All the reductions have been performed in the IRAF0 
environment, using the DAOPHOT-II routines (Stetson 
1987, 1992) in a standard way. Stars have been found in the 
deepest V frame, then every other frame has been aligned 



t Digitalization operated at STScI, on material from the UK 
Schmidt telescope, operated by the Royal Observatory Edinburgh 
and the Anglo-Australian Observatory 



t IRAF is distributed by National Optical Astronomy Observa- 
tories, which is operated by the Association of Universities for 
Research in Astronomy Inc., under contract to the National Sci- 
ence Foundation. 
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Table 3. Completeness of our measurements for the B,V,R,I fil- 
ters: this is the average of 15 experiments. Smaller prints indicate 
that the errors were higher than 0.1 mag. 



mag % B % V % R % I 



12 


100 


100 


100 


100 


13 


100 


100 


100 


100 


14 


100 


100 


100 


100 


15 


100 


100 


100 


100 


16 


100 


100 


100 


100 


17 


100 


100 


100 


100 


18 


100 


100 


100 


91 


19 


100 


100 


92 


76 


20 


100 


89 


71 


63 


20.5 


92 


76 


32 


33 


21 


80 


56 


13 


14 


21.5 


68 


28 


3 


5 


22 


50 


9 


2 


1 


22.5 


.31 


i 






23 


13 








23.5 


2 








24 


1 









Figure 2. Calibration relations: 
ous photometric bands. 



residuals of the fits for the vari- 



Table 2. Pixel and equatorial coordinates for the four reference 
stars indicated in Figure 1. 



n 


Star N. 


X,Y (pixel) 


a, 5 (2000.0) 


1 


2008 


817.16, 411.41 


16 58 49.90, -52 42 59.05 


2 


3595 


429.94, 722.28 


16 59 05.75, -52 41 05.35 


3 


890 


133.36, 191.01 


16 59 18.04, -52 44 23.63 


4 


2726 


358.67, 552.43 


16 59 08.58, -52 42 08.05 



to it, and a single list of objects was used to derive instru- 
mental magnitudes using PSF fitting. We then applied an 
aperture correction to the instrumental magnitudes, empir- 
ically derived from 10-20 isolated stars. 

Standard stars fields were analysed using aperture 
photometry. The calibration equations were derived us- 
ing the extinction coefficients for La Silla taken from the 
database maintained by J. Burki (Geneva Obs.) on t he 



U = u + 0.104(±0.033) ■ (u - i 
B = b + 0.194(±0.024) ■ {b ~ v] 
V = v + 0.016(±0.008) ■(&-« 
R = r + 0.051(±0.016) ■ (v - r 
l = i — 0.021(±0.013) • (v -i) 



- 5.219(±0.012) 
-4.640(±0.010) 

- 4.371(±0.007) 

- 4.476(±0.009) 
2.701(±0.008) 



ESO/L a Silla a rchive, accessible through www (http://arch- 
http.hq.eso.org), and are the following: 



where u, b, v, r, i are instrumental magnitudes, while U, B, 
V, R, I are the corresponding Johnson - Cousins magnitudes. 

In the central field a total of 5107 objects were detected, 
but we retained only those with internal photometric error 
given by DAOPHOT a < 0.1 mag in all bands. We then 
calibrated each frame/filter and later averaged over the same 
filter to obtain the mean magnitude in each band. For the 
shorter exposures, only the brightest bins were considered in 
the average process. Furthermore we visually inspected and 
corrected all cases where the to-be-averaged values differed 
by more than 0.15 mag: in almost all cases this happened 
for very faint objects. 

As final result, about 4200 stars have magnitudes in 
all the four B,V,R,I bands. The U exposures were much 
shallower and only about 1450 stars have been detected in 
all the 5 filters. A table with the 5 magnitudes and positions 
in pixel for all objects is available electronically from the first 
author. 

To estimate the completeness degree of our measure- 
ments we added artificial stars at random positions (using 
the DAOPHOT routine Addstar) to the deepest frame in 
each filter; we added about 10% of stars in each magnitude 
bin, distributed in colours as the real ones. The frame was 
then reduced again, using exactly the same procedure and 
the same parameters and PSF as before, and we counted the 
recovered artificial stars for each bin. The experiment was 
repeated 15 times per filter and in Table 3 we present the 
average results, except those for the U band, where com- 
pleteness reflects the much less severe crowding conditions 
with comparison to the other filters and is therefore mis- 
leading. 
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Figure 3. CMDs of NGC 6253: a) U vs U-B; b) V vs B-V; c) V vs V-R; V vs V-I. Only stars with internal errors < 0.1 mag were 
considered in the average. 



3 THE COLOUR MAGNITUDE DIAGRAMS 

In Figure | and § we present the CMDs for NGC 6253 
obtained from our reductions of the central field: in Fig- 
ure |^a,b,c,d we show the "classical" CMDs using the 5 fil- 
ters, while Figure ^,b,c is dedicated to the U band and the 
red niters. Even if the field stars dominate in number, the 
cluster main sequence (MS) is easily located in all diagrams, 
and is well defined for at least 5 magnitudes fainter than the 
Turn-Off (TO). The upper MS terminates with a "hook" 
and it shows a clear gap just below it. The TO (the bluest 
point of the hook) is located at V~14.7, B-V~0.75. The 
TO morphology is probably complicated by the presence of 
a sequence of binary stars (see later) which converges there 



in the single-stars MS. Well defined subgiant and red giant 
branches are visible, as is a red clump at V=12.7, which we 
interpret as the locus of core-He burning stars. The magni- 
tude distance between TO and red clump is 5V=2.0. Also, 
there is a well populated sample of blue stragglers. 

NGC 6325 has metallicity at the very least solar (see 
Section 4), and metal rich globular clusters (e.g. Ortolani, 
Bica & Barbuy 1993) often show anomalously extended red 
giant branches due to large blanketing effects. This phe- 
nomenon has been attributed to some open clusters as well, 
like NGC 6791 (Garnavich et al. 1994) which is very old and 
slightly more metal rich than the sun, and Cr 261 (Mazur 
et al. 1995) which is very old but with a metallicity proba- 
bly lower than solar (Friel et al. 1995, Gozzoli et al. 1996). 
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Figure 4. Further CMDs of NGC 6253 involving the U band: a) U vs U-V; b) U vs U-R; c) U vs U-I. Note the presence of a binary 
sequence above the cluster MS. 



Despite its larger metallicity, our data on NGC 6253 do not 
show such extension and turn down of the RGB, unless one 
pretends to define it through the two redder bright stars of 
the CMDs in Figure ^. We do not consider this lack as sur- 
prising, though, due to the relatively small number of red 
giants. 

From the structure of the cluster CMD we can already 
infer that it is old, but not older than some Gyr, because of 
the presence of the MS gap and hooked termination. We will 
find (see next Section) an age of about 3 Gyr. No comparison 
with previous results is possible, since this is the first study 
of the cluster. 



3.1 The Lower Main Sequence 

Note that the cluster MS looks increasingly less populated 
going towards fainter luminosities: there seems to be a flat- 
tening of the present mass function for low-mass stars. This 
of course reflects on the luminosity function and the compar- 
ison with synthetic diagrams (see next Section). This flat- 
tening is not unusual or unexpected in old open clusters (see 
the review by Friel 1995) and has already been noticed e.g., 
in NGC 752 and M 67 (Francic 1989, Fan et al. 1996), NGC 
2506 (Scalo 1986), King 2 (Aparicio et al. 1990), and as an 
extreme case in NGC 3680 (Nordstrom et al. 1995). Also 
in NGC 2243 (Bonifazi et al. 1990), a cluster very similar 
to NGC 6253, this fading of the MS is noticeable, even if 
not explicitly analysed by the authors. A likely explanation 
in old clusters is dynamical relaxation with ensuing mass 
segregation of higher mass stars towards the central parts 
and evaporation of lower mass stars (also favoured by tidal 
stripping due to encounters with massive clouds). This ef- 
fect has been noticed also in younger clusters, like NGC 225 
(age 10 8 yr, Lattanzi et al. 1991), which has photometry, 



Figure 5. Fore/background contamination as shown by the V vs 
B-V CMD of the control field. 



astrometry and spectroscopy of good quality, and even in 
very young clusters (age - 10 6 yr, e.g. NGC 2362, Wilner & 
Lada 1991), where relaxation should not have had time yet 
to work out its effects. 

Recently a similar fading of the MS has been found 
also for a globular cluster: Paresce, De Marchi & Romaniello 
(1995) presented HST data for NGC 6397, where the lower 
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Figure 6. CMDs of NGC 6253 at different radii from the centre: the cluster features, and most conspicuously the subgiant and red giant 
branches and the red clump, are well defined in all three panels 



MS, close to the H burning limit, is much less populated than 
expected. In this case evaporation of low mass stars from the 
cluster is the favoured explanation, since the relaxation time 
for the cluster is small compared to its age, and the system 
is close to the Galactic plane, even if a peculiar IMF is not 
excluded by the authors. 



3.2 Field decontamination 

The true cluster diameter is larger than that found in the 
literature (5 arcmin, Lang 1992), and one has to move 8 
arcmin from the cluster centre to find a legitimate exter- 
nal field. Figure ^ presents the CMD for the external field 
stars, extracted from the Northern frame. Given the expo- 
sure times used, we did not saturate objects in this external 
field, so the lack of stars brighter than about V=14.7 is 
real. The field comes from a mixed population, since it ex- 
hibits both a blue main sequence extending to about V=15, 
and an older/farther main sequence with TO at V=19, B- 
V=0.9 and relative subgiant /red giant branch. Given the 
galactic location of NGC 6253 and that resulting for this 
old component by extrapolating on the cluster line of sight, 
we interpret it as the old population across and beyond the 
Sagittarius spiral arm. 

Fore/background contamination is severe, but not ir- 
recoverable since cluster and field stars can be well distin- 
guished in the CMDs involving the bluer bands: the sepa- 
ration is especially noticeable in the U vs U-V, U-R, U- 
I planes (see Figure ^). Unfortunately, the U band is also 
that where less stars were detected, so we decided to at- 
tempt a cluster-field segregation in the classical V vs B-V 
diagram; there too the separation is apparent, and we think 
we can isolate bona fide field stars with reasonable confi- 
dence, at least on the blue side of the MS. Results of our 



classification are shown in Figure |8j, where the heavier sym- 
bols indicate probable cluster members, and the lighter ones 
probable field stars. All the red giants have been assigned 
to the cluster, due to their absence in the external field. 

Besides, the cluster features are well distinguishable in 
all panels of Figure [], where we show the run of the V vs 
B-V CMD with different radial distances from the cluster 
centre (taken by visual inspection of the map to be at pixel 
400,500). The number of field stars appears to decrease from 
panel a) to c) more dramatically than that of the cluster 
members, while the subgiant and red giant branches and 
even the red clump remain well delineated: we interpret this 
as further confirmation of their being cluster members. 



3.3 Binary stars 

Presence of binary stars in open clusters is not an uncommon 
phenomenon: most open clusters show in fact indications of 
a sizeable binary population. Mermilliod & Mayor (1989, 
1990) during a Coravel radial velocity survey of 10 clusters 
as old or older than the Hyades, found a percentage of spec- 
troscopic binaries around 30 % in all of them. As examples 
of photometric detections, Aparicio et al. (1990), from the 
scatter of the main sequence, found that about 50% of the 
stars in King 2 are binaries. Very well denned sequences are 
visible in NGC 2243, where at least 30% of the members are 
binary systems (Bonifazi et al. 1990) , in M 67 (Montgomery 
et al. 1993 estimate > 38% of binaries, and Fan et al. 1996 
reach 50%) or in NGC 2420 (Anthony- Twarog et al. 1990 
derive about 50% of binaries) . 

The case of NGC 6253 is not as well denned as these; we 
can clearly see (Figures ^, ^ and ^|) a secondary sequence of 
stars right above the cluster MS, at a distance of about 0.7 
mag, just what we would expect in presence of a significant 
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Figure 7. Top panels: Histograms of the colour distribution in 
U— R relative to the MS ridge line for magnitudes in the range 
U=17-20. The MS peak is at by definition, and the position of 
the secondary peak expected in each magnitude bin in presence 
of binary systems is indicated by an arrow. Bottom panel: His- 
tograms of the U-magnitude difference of all stars relative to the 
MS ridge line. The arrow indicates the expected position of equal 
mass binaries. 



appears almost exactly at the colour (indicated by an ar- 
row) where the peak of equal mass binaries is predicted to 
occur. Notice that this does not mean that the binaries have 
equal mass components: as discussed in the pioneering work 
by Maeder (1974) and most recently by Fan et al. (1996), 
binary components with mass ratios q lower than 0.5 are 
indiscernible from the single stars MS, because photometric 
errors are larger than their colour and magnitude difference, 
whereas binaries with 0.5 ^ q ^ 0.9 spread out between the 
single star MS and the equal mass binary ridge where the 
remaining binaries concentrate. In fact, both the appearance 
of the CMD and the histograms derived for NGC 6253, look 
very much the same as those shown by Fan et al. (1996) 
for M 67, where the fraction of binaries is estimated around 
50% and their mass ratios randomly distributed between 
and 1. 

An alternative test is to look at the distributions of 
magnitude differences between stars and the MS ridge line 
at each colour, since it should also present a secondary peak 
due to binaries (see e.g. Montgomery et al. 1993 for M 67) . 
The results for our cluster are not as conclusive as theirs, 
given the larger field contamination. We have computed the 
magnitude differences in the U, U-R CMD, and their his- 
togram is shown in the lower panel of Figure ^, with an 
arrow indicating the predicted location of objects 0.75 mag 
brighter than the MS: there is indeed a hint of secondary 
peak in that position, but definitely not significant. 

As it is, we cannot derive a precise figure for the bi- 
nary fraction, since we are plagued by field contamination. 
Taking into account that the number of stars in the sec- 
ondary peak of Figure histograms ranges between 0.14 
and 0.27 times the total number of stars falling in both the 
primary and the secondary peaks, the apparent fraction of 
binary systems in NGC 6253 is approximately 20%. Due to 
the various effects discussed by Maeder (1974) and Fan et 
al. (1996), this presumably corresponds to a larger actual 
percentage of binaries. 

An alternative explanation suggested by some authors 
for the MS skewness is differential extinction (due e.g. to a 
feeble dust lane, see Lattanzi et al. 1991 for NGC 225). In 
the case of NGC 6253, we reject this hypothesis, because 
the stars on the right of the MS have galactic locations uni- 
formly distributed within the cluster, and in these conditions 
differential reddening should not create a second sequence 
clearly separated from the first one as in Figure Q but rather 
a general spread of the whole MS. 



fraction of binary systems. The evidence is more convinc- 
ing when looking at CMDs built using the U and red niters 
(Figure |^), since the separation between the MS and field 
stars is larger, and the binary sequence is better discernible. 
Unfortunately, the field population smears out the MS and 
binary sequence distributions. In order to somewhat quan- 
tify the visual impression of a secondary sequence, we have 
built the histograms of the colour difference (e.g. in U-R, 
see top panels in Figure between each star and the MS 
ridge line (interpolated by eye) at different magnitude levels. 
In all panels, except the very faintest where any feature is 
wiped out by the errors, a small secondary peak to the right 
of the MS one is visible. Despite its small size, this peak pro- 
vides strong support to the binary stars hypothesis, since it 



4 CLUSTER PARAMETERS 

We have applied to NGC 6253 the approach amply described 
by Tosi et al. (1991), and already employed for two other 
old open clusters (Bonifazi et al. 1990; Gozzoli et al. 1996). 
Briefly, this method represents an improvement of the clas- 
sical isochrone fitting, allowing the simultaneous derivation 
of age, reddening and distance modulus, with the advan- 
tage of checking the relative number of stars in the vari- 
ous evolutionary phases. The observed CMDs of the cluster 
are compared with synthetic CMDs resulting from Monte 
Carlo simulations of a system containing the same number 
of stars above the same limiting magnitude, and with the 
same photometric errors and incompleteness factors in each 
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Figure 9. Synthetic CMDs derived from the FRANEC stellar evolutionary tracks. Panels (a) and (b) adopt Z=0.01, r=2.5 Gyr, 
(m-M) =10.9 and E(B-V)=0.43; panels (c) and (d) Z=0.02, r=3 Gyr, (m-M) =10.8 and E(B-V)=0.33. 



Figure 8. CMD of NGC 6253 for the objects with a < 0.1. 
Probable cluster members are indicated by filled circles, probable 
field stars by dots (see text for details). 



magnitude bin as observed in the actual cluster. For each 
assumed evolution model, the resulting synthetic CMD is 
translated into the empirical one by finding the values of red- 
dening and distance modulus providing the best agreement 
with the observed stellar distribution, and thus providing 
in turn metallicity, age, reddening and distance modulus at 
the same time. The resulting set of values, however, is not 
unique, because different stellar evolution models may lead 
to quite different solutions, as shown for instance by Gozzoli 
et al. (1996) for Cr 261. For this reason, we have derived the 
cluster parameters with three different data bases of stellar 



evolutionary tracks, and evaluated both the best parameter 
values and the corresponding theoretical uncertainties. 

The synthetic diagrams simulated for NGC 6253 are 
based on homogeneous sets of stellar evolution models com- 
puted for several initial metallicities: a) the tracks with clas- 
sical mixing length treatment of the convective zones com- 
puted by the Frascati-Teramo group (hereinafter FRANEC) , 
b) the tracks with overshooting from convective cores com- 
puted by the Geneva group (hereinafter GENEVA), c) 
the tracks with overshooting from convective cores by the 
Padova group (hereinafter BBC). We have not used the 
tracks by D'Antona et al. (1992, hereinafter CM), applied 
to Cr 261, because they are available only for stellar masses 
smaller than those describing the upper CMD of NGC 6253. 

Within the framework of each group of stellar models, 
we have performed several simulations for any reasonable 
combination of age, reddening and distance modulus, all of 
which have been compared with the empirical CMD and lu- 
minosity functions of NGC 6253. We describe below only the 
most significant cases, selected on the basis of these compar- 
isons. 

Since the synthetic CMDs obviously apply only to clus- 
ter members, we have selected from the stars detected in the 
field of NGC 6253 and with photometric error smaller than 
0.1 mag in all bands but U, only those that can be taken as 
probable members (see section 3.1). This selection restricts 
the sample to 641 objects out of 4200; all the others are pre- 
sumably back/foreground stars. The supposed members are 
represented in Figure ^ with the larger symbols, the other 
objects with the small ones. The synthetic diagrams dis- 
cussed below therefore assume the cluster to be populated 
by 641 stars. 
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We have shown in the previous section that this clus- 
ter is probably populated by a significant fraction of binary 
systems. To further check this point, we have performed our 
Monte Carlo simulations either assuming that all the clus- 
ter members are single stars or that a varying fraction of 
them are binaries. In the latter cases, for each binary sys- 
tem we have assumed a secondary /primary mass ratio ran- 
domly extracted from a flat distribution. Colour and mag- 
nitude of each system are assigned according to the mass 
ratio and following Maeder's (1974) prescriptions. This im- 
plies, for instance, that a system formed of equal mass stars 
has the same B-V as each of the two stars and a magnitude 
corresponding to double brightness (i.e. 0.75 mag brighter). 
A system with mass ratio 0.8 shows a B-V 0.04 mag redder, 
and a V 0.35 mag brighter, than the primary star; systems 
with q <0.6 have in practice the same colour and magnitude 
of the primary star. 



4.1 Results with FRANEC stellar models 



FRANEC sets of models follow the evolution of stars be- 
tween 0.6 and 1 M© in the central hydrogen and helium 
burning phases and the evolution of stars between 1 and 9 
Mq up to the onset of thermal pulses on the asymptotic gi- 
ant branch. Of the available sets with different initial helium 
and metal abundances, we have used for NGC 6253 those 
with Y and Z equal to (0.27, 0.01), (0.27, 0.02), (Castellani 
et al. 1993). | 

Figure 3 (c) shows one of the best synthetic diagrams 
obtained with Z=0.02 for a full population of single stars. 
The adopted parameters are: t~3 Gyr, E(B-V)=0.33, (m- 
M) =10.8. The shapes and locations of the MS and TO pre- 
dicted by the single star models are correct, including the 
hook, the gap and the relative number of stars in the var- 
ious phases. However, all the stars on the red side of the 
MS are missing. This inconsistency can be easily overcome 
by including a proper fraction of binary systems in our sim- 
ulations. In fact, panel (d) shows the same CMD as panel 
(c), but for a population with 30% binaries and random q, 
and reproduces quite well all the observed features of the MS 
and TO regions. We have performed similar simulations with 
other percentages of binary systems and found that lower 
fractions do not fill the MS enough, whereas higher percent- 
ages bring too many objects above the subgiant branch of 
single stars. 

A striking feature of all the performed simulations, for 
this as well as for all the other sets of stellar models, is that, 
if we adopt the incompleteness factors of Table 3, empirically 
derived from the CCD frames, the resulting synthetic CMDs 
systematically overestimate the number of faint stars. This 
phenomenon is apparent in the luminosity functions (LF) of 
Figure |l^, where the dots correspond to the observational 
data and the thick dashed curve to a model adopting the 



3 Note that these tracks have been computed with the LAOL 
Los Alamos opacities. According to the authors, the effect of 
using instead the most recent OPAL Livermore opacities cor- 
responds only to assuming a slightly larger metallicity (Cassisi 
et al. 1993). Therefore, the metallicity of the FRANEC models 
mentioned above should actually be taken as about half of their 
nominal Z value. 



Figure 10. Luminosity function for the observed (dots) and syn- 
thetic (lines) stars members of NGC 6253. Top panel: LFs cor- 
responding to the FRANEC models. The thick line results from 
adopting the actual incompleteness factors listed in Table 1, thin 
lines from that modified at the faint end (see text). The dotted 
line represents the Z=0.01 case of Figure 9 (b), the dashed line 
the Z=0.02 case of Figure 9 (d). Central panel: GENEVA models. 
The dotted line represents the Z=0.008 case of Figure 11 (b), the 
dashed line the Z=0.02 case of Figure 11 (d), and the solid line 
the Z=0.04 case of Figure 11 (f). Bottom panel: BBC models. 
The dashed line represents the Z=0.02 case of Figure 12 (b) and 
the solid line the Z=0.05 case of Figure 12 (d). 



same parameters as above and the empirical incompleteness. 
In order to predict the right number of stars in the fainter 
bins (and consequently in the bright ones as well), one has 
to artificially alter the actual incompleteness factors below 
V=17. This discrepancy is not due to an underestimate of 
the incompleteness (we have performed the artificial star 
test 15 times obtaining always very similar results), but to 
a real decrease of the number of low mass stars, beyond any 
reasonable IMF. We are inclined to interpret this decrease in 
terms of evaporation from the cluster of the smallest stars, 
similar to what Fan et al. (1996) have found for M 67. In 
all the following we then show synthetic CMDs and LFs 
(thin curves) based on the modified incompleteness which 
allows to reproduce the faint end of the observational LF. 
The adopted incompleteness start to deviate from the em- 
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pirical one below V=17, and becomes a factor of 5 more 
severe at V=21. 

Finally, the model shown in Figures [] (c) and (d) shows 
a subgiant branch more extended in colour than the actual 
branch of NGC 6253 and, as a consequence, a too red RGB. 
Binary stars obviously do not help in the RGB region, al- 
though the observed spread in the CMD distribution of the 
red giants (which, at this bright magnitudes, cannot be at- 
tributed to photometric errors) is probably due the binary 
stars effect. Due to the increased size of post-MS stars, such 
effect after the exhaustion of the central hydrogen is much 
more complicated than the simple shift in magnitude and 
colours assumed here; however, it is apparent that it cannot 
push the predicted distribution bluewards, as would be nec- 
essary to fit the data. A lower E(B-V) cannot be invoked to 
solve this discrepancy, because it would obviously improve 
the fitting of the red CMD region but worsen the blue re- 
gion. We will see in the following that this discrepancy is 
probably a metallicity effect. 

Figures ^ (a) and (b) show the best synthetic dia- 
grams, with and without binaries, obtained with Z=0.01. 
The adopted parameters are: t—2.5 Gyr, E(B-V)=0.43, (m- 
M)o=10.9. As in the previous case, a 30% fraction of binary 
stars is required to properly reproduce the colour and mag- 
nitude extension of the MS, as well as the global LF (see dot- 
ted line in Figure |ic| ) . The synthetic subgiant branch shows 
a colour extension even larger than in the solar metallicity 
case and the RGB is too red. 

To verify whether a larger metallicity could help to re- 
produce the observed cluster features, we have made some 
further simulations with the new tracks with Y=0.34 and 
Z=0.04 kindly made available by S. Cassisi and computed 
with the new OPAL opacities and a slightly different version 
of the FRANEC code (Bono et al. 1996 in preparation). This 
set is however not homogeneous to the other two, due to the 
different input physics, and this does not allow to use it 
fruitfully for our purposes. 

4.2 Results with GENEVA stellar models 

The stellar evolution tracks computed by the Geneva group 
take into account the possible overshooting of convective re- 
gions out of the edges defined by the classical mixing length 
theories. The models have been computed for the mass range 
0.8 and 120 Mq and several initial metallicities. For NGC 
6253 we have tested the cases with initial Y and Z (0.26, 
0.008), (0.30, 0.02) and (0.34, 0.04) presented by Schaerer et 
al. (1993b), Schaller et al. (1992) and Schaerer et al. (1993a), 
respectively. The low mass stellar models with solar metal- 
licity are followed up to the early asymptotic giant branch 
phase (Charbonnel et al. 1996), the others only to the tip of 
the RGB, therefore only for the solar composition cases can 
our synthetic diagrams show the red clump corresponding 
to the central helium burning of low mass stars. The lack of 
stellar models for masses below 0.8 Mq and for the clump 
makes the synthetic CMDs not completely comparable with 
that of NGC 6253, but the cluster parameters can be derived 
anyway. 

Stellar tracks with lower initial metallicity have bluer 
MS and more extended subgiant branch. The first phe- 
nomenon implies that the corresponding synthetic diagrams 
require larger reddenings to fit the observed distribution. 



Combined with the larger colour extension of the subgiant 
branch, this however leads to quite red RGBs, so that, for all 
the performed simulations, the lower the initial metallicity, 
the redder the final RGB. As a consequence, models with 
Z=0.008 and Z=0.02 are not able to properly reproduce the 
empirical CMD of NCG 6253, whereas models with Z=0.04 
are in agreement with the location of both the blue and the 
red stars. 

Figure |ll] (a) shows the best synthetic diagram ob- 
tained with Z=0.008 for a full population of single stars. 
The adopted parameters are: r=2.5 Gyr, E(B-V)=0.46, (m- 
M) = 10.9. The shapes and locations in the CMD of the MS 
and TO predicted by the single star models are correct, but 
all the stars on the red side of the MS are missing. Panel (b) 
shows the corresponding CMD for a population with 30% 
binaries. It is apparent from the shown CMDs that with 
this set of tracks, as well as with the other GENEVA sets, 
the inclusion of binary systems places too many bright stars 
right after the TO, due to the fact that the timescale of this 
phase is longer in the GENEVA than in the FRANEC and 
BBC models. We have checked whether we could reduce the 
number of binary post-TO stars by reducing the fraction 
of multiple systems, but we have not reached satisfactory 
solutions. In fact, any decrease of the binary fraction has 
much more effect on the number of MS objects than on that 
of subgiants and we end up with a clear underestimate of 
the MS binaries while still having too many early subgiant 
binaries. 

The red clump is absent because these stellar tracks do 
not reach that phase, and this absence is the reason for the 
slight underestimate of the number of bright stars app arent 
in the luminosity function (central panel in Figure |l0| ). The 
lack of stellar tracks for M<0.8 Mq is the cause the abrupt 
fall of these theoretical luminosity functions. 

Figure [ll] (c) and (d) show the best synthetic dia- 
gram obtained with Z=0.02 for a full population of single 
stars and for a population with 30% binaries, respectively. 
The adopted parameters are: r=3 Gyr, E(B-V)=0.36, (m- 
M)o=10.9. For this set of tracks the core He-burning phase 
has been computed and the corresponding red clump falls 
in the synthetic CMD at the same brightness (V~12.7) as 
the observational one, but at redder colour (A(B-V)~0.05). 
As for the Z=0.008 case, a lower reddening would not solve 
the problem, which is clearly a consequence of the intrin- 
sic shape of the subgiant branch in these stellar models. As 
for the Z=0.008 case, 30% binaries provide too many bright 
subgiants. 

Finally, Figure [ll] (e) and (f) show the best synthetic 
diagram obtained with Z=0.04 adopting: r~3 Gyr, E(B- 
V)=0.32, (m-M) = 10.9. Again, as for Z=0.008, the red 
clump is absent from these stellar tracks. Aside from this 
point, both the synthetic CMD and LF reproduce quite well 
the various data features. The Geneva models therefore in- 
dicate that NGC 6253 has a metallicity double than solar. 

4.3 Results with BBC stellar models 

Also the stellar evolution tracks computed by the Padova 
group take overshooting from convective cores into account, 
although with a treatment different from that adopted by 
the Geneva group. They have been computed for masses 
between 0.5 and 120 Mq and for several initial metallicities. 
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Figure 11. Synthetic CMDs derived from the GENEVA stellar evolutionary tracks. Panels (a) and (b) adopt Z=0.008, t=2.5 Gyr, 
(m-M) =10.9 and E(B-V)=0.46; panels (c) and (d) Z=0.02, r=3 Gyr, (m-M) =10.9 and E(B-V)=0.36; panels (e) and (f) Z=0.04, 
r=3 Gyr, (m-M) =10.9 and E(B-V)=0.32. 



Figure 12. Synthetic CMDs derived from the BBC stellar evolutionary tracks. Panels (a) and (b) adopt Z=0.02, t=4 Gyr, (m-M)o=10.7 
and E(B-V)=0.30; panels (c) and (d) Z=0.05, r=3 Gyr, (m-M) =11.0 and E(B-V)=0.23. 
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They reach the tip of the asymptotic giant branch or the 
ignition of the core C-O burning phase, depending on the 
initial stellar mass. For NGC 6253, we have tested the sets 
of tracks with Y and Z equal to (0.28, 0.008) by Alongi et 
al. (1993), (0.28, 0.02) by Bressan et al. (1994) and (0.352, 
0.05) by Fagotto et al. (1994), available at the Strasbourg 
Data center. 

As already found for the Geneva models, tracks with 
lower metallicity show bluer MS, larger colour extension of 
the subgiant branch, redder RGB. Therefore, models with 
Z=0.008 do not reproduce well the data and we do not show 
them here. Models with solar metallicity are not satisfac- 
tory either. Figures |l^ (a) and (b) show the BBC synthetic 
CMDs with Z=0.02 in better agreement with the data, re- 
spectively without or with 30% binaries. They assume r=4 
Gyr, E(B-V)=0.3 and (m-M) = 10.7. With the binaries in- 
cluded the simulation reproduces pretty well the MS and 
the TO stellar distribution both in the CMD and in the LF 
(Figure |l(j). The synthetic CMD shows too red colours for 
the giants, thus suggesting that NGC 6253 has a metallicity 
larger than solar. 

Figure |l| (c) shows the BBC synthetic CMD in bet- 
ter agreement with the data. Figure [13 (d) shows the corre- 
sponding CMD with 30% of binary stars. They both assume 
large initial metallicity Z=0.05, t=3 Gyr, E(B-V)=0.23 
and (m-M)o = 11.0 and reproduce pretty well the MS and 
post-MS colours, luminosities and stellar distributions (see 
also the LF for the case with binaries in the bottom panel of 
Figure [j^). With the Padova models we thus find, as with 
the Geneva ones, that NGC 6253 is metal rich and ~3 Gyr 
old. 



5 CONCLUSIONS 

Open clusters are a class well suited to study the structural 
and evolutionary changes in our Galaxy. In spite of all the 
uncertainties in dating them, and of the fact that even the 
age ranking is not totally safe when simply taken from litera- 
ture values, they provide ages more accurate than any other 
disk objects. Ideally one would like to determine ages con- 
sistently, fitting uniform sets of theoretical models to clus- 
ters whose metallicities and reddenings are well constrained 
by observations. This is a tough job, given the difficulties 
in obtaining precise enough photometry, and can be done 
only for a limited sample of clusters. Also, reddenings are 
usually very badly, if at all, known. Metallicities are in a 
better situation, since they have been measured for many 
clusters through narrow-band photometry or spectroscopy, 
both with low and high resolution (Friel 1995). Our method, 
with the advantage of determining at the same time age, 
distance, reddening and metallicity, can circumvent most of 
these problems and the main conclusions that can be derived 
are the following. 

a) In the case of NGC 6253 we find that age and distance 
modulus are very tightly derived with the synthetic CMD 
method. In fact, the best reproduction of the observed clus- 
ter features is obtained assuming r=3.0±0.5 Gyr and (m- 
M)o=10.9±0.1, whatever the adopted class of stellar models; 
a striking agreement if we consider that different tracks can 
lead to rather different ages, as we found e.g. for NGC 2243 
(Bonifazi et al. 1990) and Cr 261 (Gozzoli et al. 1996). 



b) The reddening derived from our analysis of NGC 
6253 is rather large, 0.23<E(B-V)<0.32, as expected from 
the galactic location of the cluster, which not only is toward 
the center, but falls on the Sagittarius spiral arm. The intrin- 
sic differences between the various sets of adopted models 
are responsible for the uncertainty on the derived value. For 
any given metallicity, in fact, the GENEVA tracks are intrin- 
sically hotter, and the BBC tracks cooler, than the others, 
thus implying systematically larger and smaller reddenings 
respectively. 

c) A good qualitative agreement is achieved also in what 
concerns the cluster metallicity. Our method cannot provide 
the precise chemical abundance of the system, since it is 
limited at least by the restricted number of chemical com- 
position cases of the available stellar tracks. However, it is 
apparent from the results presented in the previous section 
that NGC 6253 must have a metallicity roughly double than 
solar, since only in this case all the predicted locations in the 
CMD of MS, subgiant, RGB and red clump are simultane- 
ously consistent with the empirical ones. 

d) The lower main sequence shows clear indication that 
a fraction of low mass stars are missing, most probably due 
to evaporation through dynamical relaxation and mass seg- 
regation, as found in many other cases. No reasonable IMF 
could justify such an effect. 

e) The cluster has a quite distinct secondary sequence, 
right above the single-star MS, attributable to binary sys- 
tems. A sizeable population of binary stars ( ^ 20 %) not at 
all unusual for open clusters, also leads to better agreement 
between the observed and the synthetic CMDs. 

With this study we increase the number of well studied 
old open clusters, since, at an age of 3.0±0.5 Gyr, NGC 
6253 falls in this small sample, whose properties like galactic 
locations, chemical abundances, ages can be used as reliable 
tracers of disk formation and evolution. 
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